Controlling two-species Mott-insulator phases in an optical lattice to form an array of 

dipolar molecules 



O 

o 

(N 
Oh! 

(N 



O 
(N 
O 



i 

c 

o 
o 



13 



M. G. Moore t and H. Sadeghpour tt 
ITAMP, Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 

(February 1, 2008) 

We consider the transfer of a two-species Bose-Einstein condensate into an optical lattice with a 
density such that that a Mott-insulator state with one atom per species per lattice site is obtained 
in the deep lattice regime. Depending on collision parameters the result could be either a 'mixed' 
or a 'separated' Mott-insulator phase. Such a 'mixed' two-species insulator could then be photo- 
associated into an array of dipolar molecules suitable for quantum computation or the formation of 
a dipolar molecular condensate. For the case of a 87 Rb- 41 K two -species BEC, however, the large 
inter-species scattering length makes obtaining the desired 'mixed' Mott insulator phase difficult. To 
overcome this difficulty we investigate the effect of varying the lattice frequency on the mean-field 
interaction and find a favorable parameter regime under which a lattice of dipolar molecules could 
be generated. 
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In a recent experiment [1] a quantum phase transition 
from a superfluid (SF) to a Mott insulator (MI) state 
was observed by varying the depth of a three-dimensional 
optical lattice [2] superimposed onto a trapped Bose- 
Einstein condensate (BEC) of 87 Rb atoms. This ex- 
periment highlights a growing trend in ultracold atomic 
physics whereby condensates are no longer the direct ob- 
ject of study, but serve mainly as well-controlled initial 
state for the preparation of more exotic highly-correlated 
many-body states. In addition to providing new in- 
sight into fundamental phenomena of condensed-matter 
physics, MI states are expected to have important ap- 
plications in Heiscnberg-limited atom-interferometry [3] 
and quantum computing [4,5]. 

Another interesting application of a Mott insulator was 
recently proposed by Jaksch and coworkers: employing 
a Mott insulator state as an intermediate stage in the 
generation of a BEC of molecules [6] . Molecules would 
be formed via stimulated Raman photoassociation of a 
Mott insulator with two atoms per lattice site, resulting 
in a molecular MI with unit filling factor which could 
then be 'melted' to form a molecular BEC. The primary 
advantage of this approach is that during photoassocia- 
tion the resulting molecules would be completely isolated 
from collisions with other atoms and/or molecules. This 
is a distinct advantage over photoassociation in the su- 
perfluid phase [7,8], as photoassociation strongly favors 
the creation of vibrationally excited molecules- due to fa- 
vorable Franck-Condon transition from the continuum- 
which tend to undergo inelastic collisions in the pres- 
ence of other atoms and/or molecules. The resulting re- 
lease of vibrational quanta leads to loss of both particles 
from the trap. By isolating each vibrationally excited 
molecule in a lattice site, it should be possible to stimu- 
late the molecules into the ground rotational-vibrational 



state before 'melting' the Mott insulator, thus eliminat- 
ing this collisional loss mechanism. 

In this Letter, we extend this idea to the formation 
of a lattice of two-species Mott insulator phase with pre- 
cisely one atom per species per lattice site. This scheme 
would require a two-species Bose-Einstein condensate as 
a starting point, as has been recently demonstrated with 
87 R and 41 K [9,10]. The two-species condensate would 
then be adiabatically loaded into a three-dimensional op- 
tical lattice, resulting in the formation of either a 'mixed' 
or 'separated' MI state, corresponding to each lattice site 
containing different or identical species, respectively, see 
Fig. 1. The formation of a mixed phase turns out to be 
difficult for the case of a 87 Rb- 41 K system, due to the 
large positive inter-species scattering length which favors 
the separated phase. We therefore concentrate on over- 
coming this difficulty by varying the lattice frequency, 
which exploits the difference in atomic resonance frequen- 
cies to control the relative densities of the two species. 

Preparing a Mott insulator state of heteronuclear 
molecules by this approach would be an elegant method 
to obtain a lattice of dipolar molecules for use as a quan- 
tum computer [11]. In addition a rich variety of quantum 
phases 
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FIG. 1. Separated- (left) and mixed- (right) two-species 
Mott insulator phases in an optical lattice. The sepa- 
rated-phase MI is randomly distributed. 
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are predicted for the ground state of a lattice of dipolar 
bosons , including a supersolid phase [12]. The lattice of 
dipolar molecules could then be 'melted' by adiabatically 
lowering the lattice potential, resulting in a condensate of 
dipolar bosons. Here the long-range dipole-dipole inter- 
actions should introduce interesting correlations effects 
in both the ground state and collective excitations of the 
BEC [13-17]. 

Others have considered the related problem of 
superfluid-insulator transitions for spinor (multicompo- 
nent) BECs [18] in which the added freedom of superpo- 
sitions of hyperfine states is included. One crucial differ- 
ence between a two-species system and a spinor system 
is the fact that each atomic species will see a different 
optical potential as a result of having different dipole mo- 
ments, resonance frequencies, and atomic masses. This 
additional freedom should result in the physics of a two- 
species system being significantly richer than that of a 
single-species spinor system. For example, due to the 
potential for different tunneling rates between neighbor- 
ing lattice sites and/or collision rates for each species, the 
superfluid-insulator transition will most likely occur at a 
different lattice intensity for each species, resulting in a 
regime where a MI of one species is embedded in a SF of 
another species. Thus the phases of a dual-species lattice 
can be divided into three categories: a dual SF, a hybrid 
SF-MI phase, and a dual MI phase. At present we will 
focus primarily on the properties of the dual MI phase, 
as it is the only phase with the potential for collisionless 
photo-association into dipolar molecules. 

The starting point for our theoretical consideration is 
a two-species atomic system in the presence of an optical 
lattice potential. The atomic density assumed to be such 
that at the center of the trap there is one atom per species 
per cubic wavelength. The Bose-Hubbard Hamiltonian 
[2] for such a system is then given by 

H = | \3n nij{fiij - 1) + g 2 2 n 2 j(h 2 j - 1) 



where di and uii are the dipole moment and resonance 
frequency of species i, respectively, ojl is the laser fre- 
quency of the optical field, and U(r) — — \E(r)\ 2 /h is the 
species-independent intensity of the optical field. As- 
suming that the sign of detuning is the same for both 
species, in which case the potential minima will coincide, 
we can expand around the potential minima according to 



Vi(r) = i 



2 (oji - LJ L ) 



(3) 



where we have dropped the constant term consistent with 
a transformation to a rotating frame which leaves the 
Hamiltonian (1) unchanged. 

Following the expansion in (3), we make a Gaussian 
approximation for the lowest energy Wannier modes of a 
lattice site for species i as 
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where we have introduced the harmonic oscillator length 
for species i along the direction f „ 



r h 2 \Lu, 
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rrii being the atomic mass of species i. With these wave- 
functions, we can compute the collision coefficients via 
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where ciy is the scattering length for collisions between 
species i and j, and fj^j = mirrij / (rrii+mj) is the effective 
mass. 



+ 2512 hijfhj] Uhij c\ t cij + fay cl t c 2j , (1) 

where Cjj is the bosonic annihilation operator for species i 
at lattice site j, n»j = cj^c^ is the corresponding number 
operator, ftijk is the rate of tunneling of species i from 
lattice site j to site k, and gij is the collision rate for 
collisions between species i and j. 

In order to compute the various interaction and tun- 
neling rates we must first consider the dependence of the 
optical lattice potential on the intrinsic properties of each 
atomic species. The optical potential seen by species i 
can be expressed as 



Vi(r) 



(Ui - LU L ) 



(2) 
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FIG. 2. Relative density parameter in an optical lattice 
for 23 Na, 41 K, and 87 Rb as a function of the lattice fre- 
quency ujl- The relative density parameter is in units of 
[m p e 2 aQ/ci? oc ] 3,/4 and the laser frequency is in units of cR^. 
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Because the harmonic oscillator length (5) for a given 
atomic species tends to zero when approaching the 
atomic resonance frequency, tuning the lattice frequency 
close to the resonance frequency of one atomic species 
can act as a very strong 'handle' with which to tune the 
relative atomic densities, and hence the collision param- 
eters of the system. The local density per atom, defined 
as l/Vi = lip A," 1 , is given by 1/Vi = p^ 2 ^ u X J 2 , 
where we have introduced the relative density parameter 



3/4 



(7) 



From the equation for 1/Vi, we see that the relative den- 
sities of the different atomic species will be indepen- 
dent of the lattice parameters u^, provided only that 
the wells are sufficiently deep that the Gaussian ap- 
proximation holds for the lowest Wannier state- gen- 
erally required for the lattice depth to be large com- 
pared to the atomic recoil energy. In Figure 2, we 
plot the dimensionless density parameter pi as a func- 
tion of lol for atomic species currently used in BEC 
experiments: 23 Na, 41 K, and 87 Rb. From this figure, 
we note that a wide range of relative densities can be 
obtained by varying the lattice frequency across the 
atomic resonance frequencies. The atomic parameters 
in this calculation are: m* = {22.99,40.96,86.91} x 
m p , u>i = {0.9720,0.7472,0.7351} x cR^, and d; = 
{0.9036,0.8227,0.8023} x ea , respectively, where m p is 
the proton mass, e the electron charge, ao the Bohr ra- 
dius, c the speed of light, and Rco the Rydberg constant. 

In analogy with previous studies of Bose-Hubbard 
Hamiltonians [2,12,18], it is safe to assume that for a deep 
enough lattice and for strictly positive scattering lengths, 
the ground state of the two-species Bose-Hubbard model 
having N atoms per species and N lattice sites will be an 
insulator state, i.e. the atoms will be localized in individ- 
ual wells. The minimum energy arrangement of the two 
atomic species then depends on the relative strengths of 
inter- and intra-species collision parameters. In princi- 
ple, there is a broad variety of possible two-species Mott 
insulator phases. If, e.g., one species has an extremely 
weak self-interaction relative to the other interaction pa- 
rameters, then it may be energetically favorable for all 
atoms of the first species to occupy a single lattice site 
with the second species uniformly distributed over the 
remaining sites. We do not consider such extreme situa- 
tions and instead focus on the more likely regime where 
the ground state contains exactly two atoms per site. 
This regime has two phases: a mixed phase where one 
atom per species occupies each site, and a separated phase 
where each site contains two atoms of the same species. 
We note that the separated phase is unique in that the 
ground state could be any quantum superposition of dif- 
ferent patterns for distributing the two species. 

The collisional energy per lattice of the MI state is 



given by E s . m = ±(ip Sim \H\ip 8tm ) |{^ 4J }=o, where N is 
the number of lattice sites. For the case of two lattice 
sites, the quantum state of the separated MI phase is 
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whereas the ground state of the mixed MI phase is 

ivo = 4444lo>- (9) 

The energy per lattice site calculated from this two-site 
state will remain valid as the lattice is scaled up to any 
even number of sites. The energy per lattice site for the 
separated-phase is E s — \{gu +322)1 whereas the energy 
of a pair of mixed-phase per lattice site is E m = g±%. The 
condition for a mixed-species ground state is 
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FIG. 3. The ratio of the energies per particle of the sepa- 
rated and mixed phases of a 87 Rb- 41 K two-species Mott insu- 
lator as a function of the lattice frequency ljl taken in units 
of cRoo ■ 

therefore E s /E m = (gu + gTi)l^g\2 > 1, which can be 
determined with the help of Eqs. (5) and (6). We see 
that this condition differs significantly from the free space 
condition for miscibility [19]. 

In Figure 3, we plot this ratio as a function of the lat- 
tice frequency Wifor the case of 87 Rb and 41 K. We use 
the scattering length data aub-m — 99 ao, clk-k = 60 ao 
and OLRb-K — 163ao, taken from the recent report on a 
two-species superfluid [10]. In figure 3, the frequency 
domain is divided into five regions. In region I, both 
species see a red detuned lattice and the ground state 
is a separated dual MI. In region II, which spans from 
u)l = 0.733663 to the 87 Rb resonance, both species also 
see a red detuned lattice but the increase in the rubid- 
ium density results in the ground state being the desired 
mixed phase. This is because the Rb-Rb mean-field en- 
ergy is proportional to p Rb , hence the energy cost of hav- 
ing two 87 Rb atoms can greatly exceed the energy cost of 
a mixed site despite the much larger a^-K- Region III 
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lies between the Rb and K resonances, and hence the two 
species see different potential minima. Here the ground 
state will be a dual MI, but with the Rb atoms localized 
at the lattice intensity minima and the K atoms localized 
at the intensity maxima. In regions IV, which spans from 
the K resonance frequency to ujl = 0.747631 and V, both 
species see a blue detuned lattice, with region IV corre- 
sponding to a mixed phase and region V to a separated 
phase. 

The mixed phase occurs only in the vicinity of an 
atomic resonance; thus, we must make sure that spon- 
taneous emission rates are sufficiently small and remain 
negligible for experimentally realizable time scales. By 
averaging the spontaneous emission rate, proportional 
to the light intensity, over the wavefunction (4), we 

find with the help of the integral J dz j ^(iSj } e ~ Z ^ = 

6 (l ± e~ fc2 ) and Eqs. (2) through (5) that the spon- 
taneous emission rate of an atom in the ground state of a 
well of a spherically symmetric lattice with depth afiuiR, 
uor ■ = Tibj\/ {rriiC 2 ) being the atomic recoil frequency, is 
given in the Gaussian approximation by 

li = \au R —^-(l±e-^), ( io) 
2 \wi - ujl\ V ' 

where Y; t is the natural linewidth of species i and the 
+ and — signs give the result for a red-detuned and 
a blue-detuned lattice, respectively. For 87 Rb, with 
Tbj, = 5MHz, we find that the average decay rate just to 
the right of the I-II boundary is 7^ w a(l + e -1 / ^) x 

7 x 10 _2 Hz, which for a = 20 gives a lifetime per atom 
of 0.38s. Similarly, on the blue-detuned side, we find 
the average decay rate for 41 K, with IV = 6MHz, just 
to the left of the IV-V boundary is given by jk ~ 
a(l — e -1 /^) x 7 x 10 _1 Hz. For a — 20 this also gives 
a lifetime per K atom of 0.35s, roughly the same as for 
Rb. The time scale for adiabatic passage from a SF to 
a MI state is typically on the order of 10ms, so working 
in either the red (II) or blue (IV) detuned regions should 
be feasible without significant spontaneous heating. We 
note that once the mixed dual-MI state is reached, the 
lattice depth could be ramped to a large enough value 
to suppress tunneling, after which the lattice frequency 
could be moved much farther away from the atomic res- 
onance while preserving the mixed dual-MI state as a 
metastable state. With this technique the mixed phase 
lifetime could be increased significantly to facilitate fur- 
ther experimentation. 

In conclusion, we have calculated the ground state con- 
figuration of a two-species optical lattice for the case 
where the number of atoms in each species equals the 
number of lattice sites and shown that varying the lat- 
tice frequency leads to a variety of MI phases. In making 
the transition to a dual-MI state from a dual-SF state, 
it is likely that each species will undergo the supcrfluid- 
insulator transition at a different lattice intensity. While 



there are no obvious difficulties with reaching the mixed 
dual-MI state in this manner, dynamical simulation of 
the Bose-Hubbard model with the exact Wannier states 
should be performed to ensure that the system can in- 
deed be driven into the true ground state, a task we plan 
to carry out in a future work. 
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Note: During the preparation of this manuscript we 
became aware of a related manuscript by Damski et al 
[20] which reaches similar conclusions but concentrates 
on different aspects of the proposed technique. 
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